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We argue that the first order folding transitions of proteins observed at physiological chemical 
conditions end in a critical point for a given temperature and chemical potential of the surrounding 
water. We investigate this critical point using a hierarchical Hamiltonian and determine its univer- 
sality class. This class differs qualitatively from those of other known models. 
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Proteins are complex macromolecular objects jy whose structure is determined by the evolutionary process that 
formed them |2j@]. In particular this means that the folding of proteins into its three-dimensional structure must 
be efficient. There have been several attempts to grasp aspects of the protein folding process Q], in enumeration of 
configurations |Bj , in description of folding pathways j(| and in discussing the influence of water on protein structure 

0- 

Presumable proteins evolved by subsequently adding properties to simpler structures which also fold. This has 
led us to propose a hierarchical description of the folding process, where each added subunit folds conditionally on 
the folding of structures above it in the hierarchy || . Formulating this in the framework of statistical mechanics we 
suggest that the folding can be parametrized in a number of variables ip\, ifN each taking the value or 1 where 1 
correspond to a correctly folded substructure ^|J9| . The hierarchy is implemented through the Hamiltonian 



H = - £q (ipi + (fi(p 2 + p-L<P2<P3 H h <y5l<y52 • • ' <Pn) ■ (1) 

The interactions between the protein and the surrounding water may be taken into account by adding to this Hamil- 
tonian a coupling parametrized through the water variables wx, W2, Wn P|- These variables couple to the unfolded 
protein degrees of freedom because they expose hydrophobic amino acids to the water. The resulting Hamiltonian is 

H = -£q ((fi + ipiifi2 + <Pl<yS2<yS3 H 1- ^1^2 • ' ' Pn) - [(1 - Vl) w l + (1 - ^1^2)^2 + ••• + (1 - ^1^2 • • • ip N )w N ] 

(2) 

where the hydrophobic effect is taken into account by letting each of the w's take a value from the set, £ m - m + sA, 
s = 0, 1, g — 1. A is the spacing of the energy levels of the water-protein interactions. 
The partition function is 



Z = 



(e-)"(I^ +1 ) 



The variable r is the ratio of statistical weights of unfolded to folded state, per variable: 

2 1 - e~9A/T W 

with fi = —£q — £ m j n being the chemical potential of the surrounding water. 

The physical meaning of this model is that the water molecules in contact with an unfolded portion of the protein 
has lower entropy than when not in contact (thus in the our model, hydrophobicity is caused by ordering of water 
and not by repulsive potentials, as is usually believed jllj). In the model one finds that a first order transition takes 
place when the parameter r switches between r < 1 and r > 1. Plotting r against T one obtains a non-monotonic 
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function which for small [i values passes r = 1 twice, corresponding to unfolding at both low and high temperature, 
as indeed seen in experiments ]ll| , ^2| . The mechanism for the transitions is the following. At high temperature the 
entropy gain of the protein chain causes the unfolding. As temperature is lowered the system gains more entropy by 
shielding the hydrophobic residues from the water. This leads to folding. As the temperature is lowered even further 
the cold unfolding transition occurs. Below this transition entropy is insignificant and the dominating effect is the 
attractive coupling between the water and the unfolded protein. 

For an intermediate value of the chemical potential, r just touches the line r — 1, that is dr/dT — when r = 1, 
corresponding to a merging of two first order transitions. This defines a critical point. Around this point, r varies 
quadratically in T — T c and linearly in /i — [i c , as seen from expanding Eq. (|J). In experiments of protein folding this 
point is accessible by changing the pH value of the solution. In fact, Privalov's data on low pH values indeed indicate 
that such a critical point exists. The scaling properties around this point thus opens for a possibility to gain insight 
into the nature of the folding process, in particular whether the hierarchical scheme we suggest can be falsified. 

In Fig. [l]a we show heat capacity as a function of temperature for chemical potential below, at and above the critical 
value = fi c . For the chosen values of £q = 1 and level density A = 0.02 and g — 350 the critical point is situated 
at T c — 1.33303 . . ., fj, c = 1.2838 . . .. That is, it is situated at a minimum of the heat capacity curve. This is at 
first sight surprising, usually heat capacity has a pronounced increase at the critical point. The minimum reflects a 
partial ordering of the hierarchy, as envisioned in Fig. [l]b where we show the degree of folding, counted by the average 
number of folded variables ipi = 1, i — 1, n from i = 1 until the first variable i — n + 1 which takes value tp n + 1 = 0. 
The average value of this (n) is N/2 at the critical point, reflecting that the system is on average half ordered at this 
point. Correspondingly the heat capacity dips to a value in between the value of an unfolded and a completely folded 
state. 

To characterize the functional form of the dip in the heat capacity, we investigate analytically C S i„ g (T) — C(T, fi) — 
C(T, fi c ) with fx >> fi c for different values of hierarchy size N. For finite N we may express the singular part of the 
heat capacity in the form: 



where g(x) — > const when x — ► oo and g{x) oc x a when x — > 0. We find analytically a = v = 2 from differentiating 
the partition function (|^). Fig. ||a demonstrate this finite size scaling. Similarly we in Fig. ^|b show the behavior of 
the order parameter (n) as function of T — T c and N: 



with f(x) — > const when x — > oo and f(x) oc x~@ when x — > where exponents (3 — —2, also found analytically. It 
may be surprising that /3 is negative, but this reflect in part the unusual use of an extensive (in N) order parameter, 
in part that for fi = fi c then the order parameter only obtains a non-zero value at T = T c when N — > oo. 

Likewise, we find that the susceptibility x = d(n) /dfi scales as \T — T c \ 1 where 7 = 4 and that (n) oc (jjl — ^c) 1 ^ for 
/i > /i c where 5 = — 1. Thus the usual exponent relations, a + 2/3 + 7 = 2, a + (3(6 + 1) = 2, and 7(5+1) = (2 — a) (5 — 1) 
are fulfilled Jl3[ ] . However the hyperscaling reation dv = 2 — a, where d is the dimensionality of the system, is not 
fulfilled. However, this relation has no meaning, as there are no spatial degrees of freedom. 

In terms of experiments on proteins, the relevant scaling behaviour is the how the degree of folding (order parameter) 
and the heat capacity behaves as function of temperature, when one changes chemical potential away from its critical 
value. The qualitative prediction is that the width of the singular part of the heat capacity has a minimum at the 
critical value \x = /i c . The broadening of the heat capacity is 



where h(x) oc x~ 2 for x — * oo and h(x) = const for x — > and where A/i = max(/x — /i c , A/i m j n ) with Afi min oc \/N 
takes into account the finite size sensitivity of the scaling. We show in Fig. Bp,, an example of such a data collapse. 
These predictions are experimentally accessible through the use of standard calorimetric techniques, where one should 
seek to obtain a data collapse above the critical point, i.e. the point of minimal width. The heat capacity below the 
critical fi is complicated by the merging of two first order transitions. However, the distance between these moves 
away from each other in T as A^ 1 / 2 . 

Likewise, we expect the degree of folding (n) to show data collapse of the form 




(5) 




(6) 




(7) 
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(n){T-T c f = k 



( 



T-T c 



) 



for [i> fi c 



(8) 



where k(x) behaves asymptotically as h. We show this in Fig. |3Jb. This quantity can be observed experimentally 
through fluoresence measurements. 

In summary, we have proposed that the hierarchical ordering of folding pathways implies a critical point with a 
diminishing heat capacity at criticality. We have determined all critical exponents, and proposed two experiments 
that could confirm or falsify the concept of hierarchical folding. 



[1] B. Alberts, D. Bray, J. Lewis, M. Raff, K. Roberts and J. D. Watson, Molecular biology of the cell (Garland Publ., New 

York, 1994). 
[2] H. Neurath, Science, 224, 350 (1984). 
[3] W. Gilbert, Science, 228, 823 (1985). 

[4] K. A. Dill, S. Bromberg, K. Yue, K. M. Fiebig, D. P. Yee, P. D. Thomas and H. S. Chan, Protein Science 4, 561 (1995). 

[5] P. A. Lindgard and H. Bohr, Phys. Rev. E 56, 4497 (1997). 

[6] K. A. Dill, K. M. Fiebig and H. S. Chan, Proc. Natl. Acad. Sci. 90 1942 (1993). 

[7] H. Li, C. Tang and N. S. Wingreen, NEC preprint (1997). 

[8] A. Hansen, M. H. Jensen, K. Sneppen and G. Zocchi, Physica A, in press (1998). 

[9] A. Hansen, M. H. Jensen, K. Sneppen and G. Zocchi, submitted to Proc. Natl. Acad. Sci. 
[10] J. T. Edsall, J. Am. Soc. 57 1506 (1935). 
[11] P. L. Privalov, Biochem. and Molecul. Biol. 25, 281 (1990). 
[12] G. I. Makhatadze and P. L. Privalov, Adv. Prot. Chem. 47, 307 (1995). 

[13] H. E. Stanley, Phase Transitions and Critical Phenomena (Cambridge Univ. Press, Cambridge, 1971). 

FIG. 1. a) Heat capacity, C, as a function of T. b) Degree of folding, (n), as a function of T. Here g — 350, A = 0.02 and 
N = 100. The value N = 100 has been chosen as to be close to realistic values for this parameter. 

FIG. 2. a) Finite size scaling of the heat capacity for [i = fi c , g — 350 and A = 0.02. Here a = 2 and v = 2. b) Finite size 
scaling of degree of folding, (n). Here j3 = —2. 

FIG. 3. a) C stng (T ~ T c ) 2 vs. (T - T c )/A^ 1/2 . b) (n){T - T c ) 2 vs. (T - T c )/A^ 1/2 . We have chosen N = 100, g = 350 and 
A = 0.02. Note the good quality of the data collapse in spite of smallness of the system. 



3 




Figure la 

A. Hansen, M.H. Jensen, K. Sneppen and G. Zocchi 
Critical aspects of hierarchical protein folding 




Figure lb 

A. Hansen, M.H. Jensen, K. Sneppen and G. Zocchi 
Critical aspects of hierarchical protein folding 




Figure 2a 

A. Hansen, M.H. Jensen, K. Sneppen and G. Zocchi 
Critical aspects of hierarchical protein folding 




Figure 2b 
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